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1 Introduction 



In the gauge-gravity duality picture the gravity dual of finite temperature = 4 supersym- 
metric Yang-Mills theory is a 5-dimensional AdS space (times S5) with a black hole PHI2]. 
While this is a very special conformally invariant theory and requires also g'^ Nc ^ I, Nc ^ 1 
for its validity, it has recently become apparent that one could meaningfully apply the gauge- 
gravity duality picture to QCD matter produced in relativistic heavy ion collisions, to vis- 
cosity O HIS], to jet energy loss [6l [3 El [9l [lOl E] and to photon production [12j. Note that 
these are all quasistatic phenomena; the gravity dual of a full dynamic heavy ion collision is 
unknown, for attempts, see [131 [HI [T5] . 

Ideally, one would like to derive the gravity dual of QCD in a string theory framework, but 
this top-down way has so far not lead to a unique result. Instead, one has, in a phenomeno- 
logical approach, studied various deformed metrics as candidate gravity duals of QCD-like 
theories [I6l[171[l8l[ll|20l|2ll|22l|23]. The purpose of this article is to continue this line 
and to present a gravity dual model for QCD thermodynamics. The model will cover both 
the deconfined phase (stable at T > Tc, metastable at T < Tc), the confined phase (stable at 
T < Tc, metastable at T > Tc) and the phase transition in between. The model for the high 
T phase will be a deformed AdS black hole metric with a dilaton, the model for the low T 
phase will be a deformed horizonless AdS metric with a dilaton. A deformation is defined as 
a multiplication of the 5d metric by an arbitrary function depending on the fifth coordinate 
z. The dilaton plays an essential role here. The phase transition will be analyzed in two in- 
dependent ways: by computing the energy momentum tensor in the 4d boundary theory and 
by computing the difference of the 5d bulk actions [241 [2l [25| [26] . The shapes and parameters 
of the deformations of the two phases are fixed by using the QCD number for the Tc of the 
4d boundary theory and by using constraints from the T = hadron spectrum [23]. With 
these parameters the 5d bulk Hawking-Page transition coincides with the 4d boundary one. 
With the same parameters one can then also compute the QQ free energy (see also [271, l28] 
in the two phases. The significant role of the dilaton is apparent here, the deformation in 
the string frame is qualitatively different from that in the Einstein frame. At T = we see 
that the QQ potential contains a confining linear part but also string breaking, indicating 
that the deformation contains effects of dynamical quarks. At T > Tc the QQ pair becomes 
deconfined with entropy ^ 2.1/pair. 

The model we end up with is admittedly not unique and contains several parameters, but 
this is due to the fact that we wish to treat a large number of different QCD phenomena at 
the same time. Limiting oneself to one phenomenon at a time, simpler gravity duals can be 
presented. Our result shows that a comprehensive solution is possible. 

The paper is organized as follows. In section 2, we specify the model. We assume a five 
dimensional gravity action with a dilaton and a matter Lagrangian and postulate two solu- 
tions corresponding to QCD-like theories on the boundary. In section 3, using gravity/gauge 
theory duality we compute the energy momentum tensor on the 4d boundary and find the 
phase structure of the boundary gauge theory. In section 4, we use the meson mass spectrum 
to constrain the parameters of our approach. In section 5, we consider a Hawking-Page type 
transition of the 5d bulk geometry and find the correspondence between the phase transi- 
tion on the boundary and the phase transition in the bulk. In section 6, we compute the 
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Wilson-Polyakov loops and correspondingly the free energy (potential at T = 0) between 
a probe quark and antiquark. This confirms that, indeed, at high temperatures we are in 
the deconfined phase and at low temperature in the confined phase. Section 7 contains a 
discussion. 

The numbers related to experimental ones quoted below should be understood to have a 
sizable 0(10%) error. 



2 The Model 

We consider the following five dimensional action in Einstein frame, in standard notation: 

« = lek / "'"^ ('^ - - ^''•*'') - s4 / '^'^^ + / <''-^»°*^- 

The second term is the boundary term (for a boundary surface z = e) and the third term is 
the matter part, which we do not have to specify completely. The constant a is dependent on 
the coupling of matter to gravity. This action can be the bosonic part of a supergravity action 
descendant from string theory. Note that even though we do not specify the origin of the 
model in detail, we have in mind near horizon geometry of stacked branes solution that gives 
asymptotically AdS^ space. A nontrivial matter part can be introduced if we add, for example, 
a spacetime filling brane adding fundamental matter to the problem j29 t 130 ^ IHTl I32j . The effect 
of the presence of fundamental matter would be already reflected in the full supergravity 
solution. 

We clSSlllXlG thclt Lrn 

can be so chosen that the following metric is a solution as a dual to 
finite temperature QCD. 



h{z) 



Z^ 



dz'^ 

-g{z)dt^ + dxl + dxl + dx| + 



, (2) 

with a dilaton ^{z). The function h{z'^), h{0) = defines the deformation. The boundary 
theory lives at z = and the geometry would be asymptotically (when z — > 0) AdS. We shall 
define the string frame so that it contains the combination ^J—ge~^^R in the action. The 
transformation between Einstein frame and string frame is then given by 

{E) _ _i# (s) 

9mn — ^ ^ 9mn l-^J 
and the bulk metric in the string frame will be the same as ([2]) but with 

eM-)^e^(^), f{z)^h{z) + U{z). (4) 

Our gravity dual is defined by specifying the deformation h{z) and the dilaton ^{z). We 
will have two alternatives, one will be the stable phase for T > Tc, the other for T < Tc- The 
former is 



h{z) = —cz^ 
3 
4' 



Mz\ = -6z^ 



9{z) = 1-4' (5) 
Zq 
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and the latter 



hiz) 
9{z) 



-cz\l 
-cl>z 

1. 



29 
— 1 
20 



-w 



(6) 



Here 



1 



(7) 



is the horizon radius of the metric ([2]), c ~ 0.127 GeV^ is a constant determined from the Tc 
of the QCD phase transition in the 4d boundary theory and cp ~ 0.311 GeV^ is a constant 
which, together with w{oo) = w{x — > oo) = —1.0 is determined from the mass spectrum 
of QCD. A further property of the function w{x) is that it for small x should behave as 
w{x) = + also constrained by the QCD phase transition in the 4d boundary theory 

(to make the bag constant vanish in the low temperature phase). We can thus constrain the 
UV and IR limits of ■w{x); the intermediate behavior affects details of the QQ potential. We 
will use the ansatz 

w{x) = 2x^6""^' - tanhx^ (8) 

plotted in Fig(TJ In the same figure we also show the combination f{z) = 29[w{^/cz) — 
w{oo)]cz'^ /20, later identified as the deformation of the low T metric in the string frame 
Eq. ()45p . Further important constraints come from the fact that we want the phase transition 
in the 4d boundary theory and the Hawking-Page transition in the 5d bulk theory to agree. 
In particular, this is possible only if the metric deformations coincide for z ^ 0, i.e., the 
constants c and (p coincide in ([5]) and The 5d bulk transition further constrains c > 

0, ■w{(X)) < 20/29. Also note in the ansatz ([5]) and we have assumed the dilaton and 
matter from Lm have the same leading order behavior. 





z [GeV 1] 



Figure 1: A plot of w{z) defined by Eq.® (left) and a plot of the low T metric string frame 
deformation f{z) = 29[w(V^z) - w{oo)]cz^ /20, c = 0.127 GeV^ u;(oo) = -1. 
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3 Thermodynamics of the Gauge Theory on the Boundary 



We know that the metric ([2]) for /i = is a solution of the 5d AdS equations 

1 6 

Rain - -^Rqmn - -^qmn = (9) 

and the boundary theory is hot = 4 supersymmetric Yang-Mihs theory in 4d at the 
temperature T = I/ttzq with a pressure = 3/4 times the ideal gas value [33] . How is this 
affected by the introduction of the deformation h{z)'! To be able to treat the two cases 
(0) and ([HD simultaneously, write h{z) = diz"^ + d2z'^ + 0{z^)^ only the parameters di^d2 of 
dimension T'^,T^ matter for this problem. 

To answer this question we should compute the energy-momentum tensor of the boundary 
theory. Since the theory is strongly coupled, the computation directly from the boundary 
field theory is not feasible, but we can use the gravity /gauge theory duality for this purpose. 
How this is to be done has been studied, say, in [Ml [351 ESI EZl [38]. However, the general 
methods involve writing down an explicit 5d gravity-matter action and obtaining the metric 
and other matter fields as solutions of equations of motion. Here our starting point is a 
phenomenological deformed metric and the gravity action contains an unspecified matter 
Lagrangian. This unknown part will lead to an unknown part in the energy momentum 
tensor we obtain. It will cancel when studying the phase transition, computing the difference 
of pressures between two phases. Analoguously, in the later (Section 5) analysis of the phase 
transition as a bulk Hawking-Page transition in 5d, the unknown matter Lagrangian can be 
eliminated using equations of motion. 

We transform the metric ([2|) to the form 

ds^ = -^(gnuix, z^dx^'dx" + dz^), (10) 
where g^i, r]^i, at z = 0. If we expand the metric gf^i^ for small z: 

g>.u{~z,x) = 7?^,(x) + g'j^J{x)z^ + gl^J{x)z^ + gj^JixYz^ + • . . , (11) 
the energy-momentum tensor T^,y of the boundary theory is given by [37] 

T,Ax) = 4^[#J(^) + + V] = + + Y,,], (12) 



where 



= -g^A*<^[(ti^5(2)) -tr5(2)]- 2 5(2)^^+4*1^ 5(2)- 5(2)^!^ 
= -C'^^i.v if g{2)iiv = Cr]^,y (13) 

depends only on the square of contractions of g{2) and where 

V = CmTii^u (14) 



4 



is a part depending on matter fields. We also used that for AdSs x spacetime 



3r5 



GlO = G^TT C 



vr 



t£8 



2A^2 



(15) 



For a thermal stationary and homogeneous system one, of course, has g^uiz,x) = g^y{z). 

Although we cannot determine the constant Cm in ()14p . we expect it to be the same for 
both phases. This is analoguous to the fact that in our metric ansatz ([5]) and ([6]) the leading 
boundary behaviour of both the deformation h{z) = —cz^ + • • • and the dilaton '^{z) = 
is the same in both metrics; which also imply that this holds for g{2)nu computed below. 

For the deformed metric ([2]) the transformation relating z and z is defined by 



dz 
z 



(16) 



where we have enforced the condition z ~ z for z ^ by putting the same e as the lower 
limit in both integrations. This can be integrated to 



2^2 



~ — exp 
1 + V^ 

For h = 0, unmodified black hole, one can invert 



dz 



oh{z)/2 _ ^ 

z\/g{z) 



(17) 



1 + zV(44) 
and get for the metric the exact result of |15j 



£2 



ds = 

z^ 



i\2 



Z^ 



-ii-^^dt^ + (1 + ^){dx\ + dxl + dxl) + dz^ 



(18) 



(19) 



For any /i, to invert (jlTh to express z^ in terms of z^, we expand in z^ using the expanded 
form h{z^^ = d\z^ + diz'^ + higher terms do not matter. One then obtains, for the high 
temperature metric ([5]) with the horizon zq 



g,iu = Vfj-u - ^dirjf.^z'^ + [{u^Uu + ^V/iu)^ 



20d2 



7d^ 



16 



(20) 



where = (1,0,0,0) so that u^Ui, + l?]^^ = diag(3/4, 1/4, 1/4, 1/4). Using (fT2|) .(fT3]l.(fT 
and (lisp the energy-momentum tensor then is, for the high temperature metric. 



where 



2 at2 



( SaT^ + Bh Q 

aT'^-Bh 

aT^-Bh 

\ aT"^- 



Bh J 



(21) 



tt'N 



_Nl_ 
32^' 



N 



-20d2 + 29di - l6Cm) = ^(29c' - 16C7^), 



f22) 
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where we have inserted for the deformation of the high temperature metric ^ di = c, d2 = 0. 
For the low temperature metric ([6]) there is no term, 

Tf,u = -Bir^^u (23) 

with 

Bi = ^(-20(i2 + 29dl - 16Cm) = ^(-'^^Cm), (24) 

where we have constrained the low temperature metric deformation by 20^2 = 29di. The 
motivation for the omission of the term is that in the QCD phase transition there is a 
large change in the number of degrees of freedom when one crosses Tc- In fact, at large 
A^c, p'h(Tc) ~ 0{N^) while p^(Tc) ~ C(l) [39]. In the high T phase there are gluons and 
quarks while the low T phase only contains massless pious. The motivation for the constraint 
20d2 = 29di is that it minimises the vacuum energy in the low temperature phase in which 
the hadrons move in the true physical vacuum. For example, in the bag model for hot QCD 
matter one writes Ph{T) = UhT^ — B, pi{T) = aiT'^ with ^» ai 0. 




Figure 2: The pressure p{T) and entropy density s{T) = p'{T) in terms of temperature T with 
B = 29c27V2/(327r2) and L = T^As = 2%Nlc'^ / {Stt'^). The contribution of the matter fields is neglected 
here, Cm ^ c^- 

Thus our result for the pressures of the two phases is 

2 Ar2 

pu{T) = —^T'-^{29^-mCm) 

PiiT) = ^(16C™), (25) 



and for the entropy density 



s,{T) = 0, (26) 



Though the lattice data rather follows the pattern p = UhT — BT [40] 
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The phase transition occurs at pi{Tc) = ph(Tc). The unknown matter field contribution 
Cm cancels from here, leading to 




This is what one dimensionally expects, introducing a deformation exp(cz^ + • • •) one has, 
in addition to T = I/ttzq, introduced a new energy scale c^/^ to the problem. Using ()27p to 
determine the constant c we have, e.g., 

rc = 186MeV^c« 0.127 GeV^. (28) 

The general behavior of the pressure p(T) and entropy density s{T) = p'{T) is shown in Fig. 
[21 The transition corresponds to first order phase transition with two metastable branches 
persisting at all T. In section 6 we will show that the stable branch T > corresponds to a 
deconfined phase and that at T < Tc to a confined phase. 

For high temperature geometry ([5]), one can get the entropy directly from the geometrical 
data of the bulk, computing the area of black hole or evaluating the Euclidean action. If we 
calculate the entropy density from the area law, we get 

Vr^iV^ n 3 C 

SH{T) = ^r^e-^^. (29) 

This looks very interesting phenomenologically since it fits well with the lattice data[40] with 
proper adjustment of number of degrees of freedom. However, it differs from the result (j26p 
for the entropy density. This is again due to the fact that we start from a deformed metric 
and do not derive it as a solution of 5d gravity- matter equations. For this complete treatment 
the results should concide [H] changing also the black hole entropy formula ()29p . 



4 Meson Mass Spectrum 

In this section, we show how it is possible within our framework to get a mass spectrum for 
mesons which grows linearly with the excitation level n and spin 5, i.e., rn? ~ in + S). We 
use the constraints, considered first in ^3]: in the limit z — > oo, the function f[z) and the 
dilaton ^{z) should behave as f{z) z^ , where k < 2, and ^{z) z"^, to obtain a linear 
mass spectrum as a function of n and S. 

At zero temperature, we have the metric ([6]). We specify the matter part in ([T]) to contain 
a spin S field : 

S = J d'x^e'" (^^{DAs? + Im'aI^ , (30) 

where A is totally symmetric rank S tensor field. We are assuming that the matter in this 
case is open string matter on the spacetime filling branes, for example, and the coupling to 
the dilaton is fixed in this spirit. In the axial gauge, the equation of motion for the transverse 
traceless part, after rescaling As = e^^~^^^^^^ As with F{z) = f{z) — 21og(2;), is 

(e(^-^)^e-*a.i5) + m2e(^-^)^e-*is = 0. (31) 
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After the substitution As = ip with B{z) = ^{z) — {S — ^)F(z), the equation is reduced 
to the Schrodinger form 

(-af + ip = (32) 

with the potential 

(33) 



V{z) = -(d,B 



)^ - \dlB. 



To have the expected behavior of mass rn? ~ n, the asymptotic form of the potential should 
be quadratic [23]. This constrains the form of the metric strictly. We can read off of the 
asymptotical mass spectrum of spin S meson from (|2l6p . 



m 



(S-\) 



29 

c(l wico)) 

^ 20 ^ 



(n + 5), 



(34) 



where w{oo) is the limiting value of the function w{z) in ([6]) when z oo. We have also 



used ^{z) 



')Z^ 



To match the experimental data [42], the term multiplying S" — 1/2 must 



be canceled. This fixes the value of the dilaton: 



c (l - ^u;(oo)) . 



(35) 



For the n-excitations of the rho meson with spin S 
this cancellation [fl 

3(/)n. 



A linear fit to the observations is 



1, IGJPC = 1+1- 

one obtains after 
(36) 

= 0.31 GeV2 



0.93n GeV^ [32] 

and, using the earlier established value c = 0.127 GeV^, 
It is natural to expect that ^y^f ~ e**^^ . In our case 
strong coupling regime from UV to IR. Even though the direct comparison of the weakly 
coupled and strongly couple regime is impossible, this is still qualitatively consistent. In our 
notation, this also gives the condition that u)(oo) < ||. 



Eq. ([36|) then gives (j 
35l) gives w(qo) = —1.0. 
> 0, and the coupling flows toward 



5 Hawking-Page Transition of the 5d Bulk Geometry 

In section [3l we have computed the boundary energy momentum tensor for the deformed 
geometry which corresponds to finite temperature QCD on the boundary in terms of grav- 
ity/gauge theory duality. Notice again that the energy momentum tensor (|2ip is that of the 
4d boundary theory and the phase transition is based on the 4d boundary theory picture. 
In this section we will try to compute the phase transition in the 5d bulk a la Hawking and 
Page [Mjfl- First we Euclideanize the metric ^ performing the Wick rotation t — > ir, then 
evaluate the Euclidean action ([T|) corresponding to the two metrics ([5]) and ([6]). The metric 
corresponding to the smaller action is the stable one. 

^The near boundary structure of the function w{z) affects the exact mass spectrum f53], and therefore 
needs a closer analysis. 

''A similar idea was used for the computation of deconfinement temperature in [26] for a deformed AdS 
geometry. 
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The action ([T]) contains the unspecified matter term Lm , but actually this can be eliminated 
using the equation of motion of the dilaton. The contribution to the action from the matter 
part is given by 

= / d^'xd^i^g'^d,^] = / dTd'x\lzTn^9''dM (37) 

The upper limits z = zq and z = oo correspond to the metrics ([5]) and ([6]) and both give 
a vanishing contribution. The contribution from z = e is divergent but the divergence is 
the same for the two metrics and thus cancels when the difference between the actions is 
evaluated. Thus the Lm-term can be entirely omitted from the discussion. 

To evaluate the action for the remaining terms we have to specify the limits of integration 
in / d^x = J drdzd^x. The integration over x is over all the space and thus gives only the 
volume factor V3. The integration over z is over (e, zq) for ([5]) and over (e,oo) for ([6]). Since 
we are assuming that the black hole is in thermal equilibrium with thermal radiation, the high 
temperature geometry ([5]) has a natural temperature defined by the black hole temperature 
Ph = t^zq and the integration over r is over the interval (0,/3/i). The temperature of thermal 
radiation in ([6]) can be arbitrary. For comparison with the high temperature metric, the 
temperature is fixed [24^ [2] requiring that both have the same physical circumference along 
the Euclidean time direction r at the regularization point z = e. This leads to the relation 



(38) 

Since the r-integrals just give the upper limits, the action difference to be evaluated now 

is 

A5^ = 5f -5f = -^^(^£"dzL;,-£^"dzLi) +A5boundary, (39) 

where Lh and Li are the Lagrangian densities without the L^-term in ([T|) evaluated with the 
metrics ^ and ([6]), respectively. The z-integrals in (f39]) converge since c > and since the 
earlier result w{oo) = —1.0 satisfies w{oo) < 20/29. The e — > divergences in ([39|) require a 
more careful analysis. The divergent parts of both of the Lagrangians are the sam^: 



Ld^v = -4 - ^ + (27c' - 302)1. (40) 

z° z-^ z 



but one also has to take into account the e term in Pi/Ph^ Eq. ([38j) . The result is: 



A5^' = ( dz{Lh - Li) - dzLi - C ) , (41) 



where C = -203cV20 + l/z^ denotes the sum of the correction terms coming from the 
boundary term (—580^/5) and the temperature matching (+290^/20 + 1/zq). 



^This can be guaranteed by choosing the same value for the parameters (j) and c, in both metrics. 
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Computing numerically the difference of the actions for the metrics ([5]) and ([6|) and using 
the parameter values obtained earlier from the 4d boundary theory and the function 'w{z) in 
([8]), we obtain the result plotted in Fig. [3l in arbitrary units. The phase transition occurs 
when AS^ = and this corresponds to Tc ~ 0.181 GeV. This is essentially the same as the 
Tc = 0.186 GeV used to constrain the parameters of the 4d boundary theory. Exact agreement 
can be obtained by somewhat modifying the function w{z) in ([8j), for example, by writing 
w{z) = 2z'^e~^'^'^^^ — tanh z'^. Again we emphasize that we use the bulk computation to match 
the critical temperature obtained from the boundary computation instead of comparing the 
full thermodynamic functions to that of the boundary theory in every detail. 

1 

0.5 



CO 

-0.5 

-1 
-1 . 5 

-2 

0.1 0.2 0.3 0.4 0.5 0.6 

T [GeV] 

Figure 3: The difference Sj^ — between the actions computed with the metrics ([5]) and The 
difference vanishes at T « 0.181 GeV, above that Sj^ < Sf" and the high T metric ([S]) dominates. 

6 QQ Potential and Confinement-Deconfinement Transition 

The QQ potential (or free energy at finite T) in QCD can be "experimentally" measured 
using numerical lattice Monte Carlo techniques. Basically, one measures expectation values 
of rectangular Wilson-Polyakov loops with one of the sides becoming large. At finite T, the 
correlator of two Wilson-Polyakov lines {Tr W (L)Ti W {0)) = exp{—F{L,T)/T) measure the 
free energy at finite T. 

To evaluate the QQ potential we use standard gauge theory /gravity duality techniques (see, 
e.g., [Ml US]) to evaluate the finite T expectation value of a Wilson-Polyakov loop of quarks of 
fundamental representation. Consider first a temporal-spatial loop. The Nambu-Goto action 
of a fundamental string then is 

Sng = ^ / drda^J det GMNdaX^'^dbX^ . (42) 

Since we are interested in a temporal Wilson loop we choose the string coordinates as = 
{t, x,0,0, z), where z = z{t,x) = z{x) is independent of time. The end points of the string 
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are fixed at 2;(0) = ibL/2 and z{—x) = z{x). With the gauge choice of r = t and a = x the 
string action becomes 

rL/2 pf{z) 



•S'ng 



£2 



27ra' 
£2 



dx- 



At 



Tra' 



-L/2 

z,(L) (,f{z) 



dz- 



■^z'{xY+g{z) 
^l+g{z)x'{z)\ 



(43) 



where z^, is the maximum value of z (Fi^U and e is a parameter used to regulate the UV 
divergence at z = 0. A potential will be defined by 5ng = V[L)/S.t. 

The string frame deformation f{z) is given in terms of the Einstein frame deformation h[z) 
by Eq.(jl]), f{z) = h{z) + 4$(z)/3. Using Eqs. dSD, ^ and ([35]) we have 



29 

f{z) = {(j) — c)z'^ = —[—w{oo)\cz'^ highTmetric, 



'^[wi^Jcz) - w{oo)]cz^ 



low T metric, 



(44) 
(45) 



where we earlier fixed w{oo) = — 1 and c = 0.127 GeV^. Due to the rapidly growing dilaton 
term the string frame deformation thus behaves qualitatively differently from the Einstein 
frame one: it grows monotonically ~ +z'^ for the high T metric and is always > for the low 
T metric. 

The equation of motion following from (j43p is simple and its first integral is 



of(z) 



z^ ^z'{xY+g{z) 



^f{z,) 

constant = — —\Jg{z^), 

z^ 



(46) 



where the constant is given as the value of the left hand side when z'{x) = 0. From the 
symmetry of the problem this is at x = 0, i.e., z^, is the maximum value of z. Basically the 
solution can be of the three types as in Fig. [H 






z = 



Figure 4: Three types of solutions: a QQ solution (left), a solution with z' = Q only a,t z — z^, 
(middle), a solution with a constant piece &t z = z^ (right). 

At high temperature with the metric ([5]), the equation of motion (j46p has the following 
types of solutions (Fi^: 



A quark-antiquark solution (Fig|4]left) 

x = ±L/2, x'{z) = 0, z'{x) 



oo, 



zo, 



(47) 
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with the action and the corresponding potential, regulated as discussed below, 

£2 



Vn 



20 ^fi^) 

dz TT- 



ira 

vra' 
£2 r 



1 



+ 



z^ 



z^ 



1 



-rra' Uo 



,/{2) 



1 



1 



(48) 
(49) 



A solution connecting the quark and antiquark with z'{x) = only at a; = (Fig|4] 
middle). Integration of (06]) relates z^, and L: 

1 



L 



dz 



where 



qiz) 



Voiz) [q{z,)/q{z) - 1] 



-2/(2). 



Zn - 



for high temperature phase and 



for low temperature phase. The corresponding potential is, regulating the z 
divergence as discussed below, 
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One solution of the equation of motion also is (FigU] right) 

z = = constant. 



^ ] z* 



(53) 
(54) 

(55) 



This can be made a solution of required type by connecting its end points to the points 
2;(0) = ibL/2 by a solution of the previous type with nonzero z'{x). If the length of the 
constant piece is £, the value of z^, will then be determined by (|50p with L ^ L — i and 
the total potential is evaluated for this z^ (L — i) added to 



£2 e/(--) 



planar " o 



\fg{z7) t 



(56) 



For fixed L = ^ + (L — £), < £ < L, we thus obtain a family of solutions with increasing 
length i of the constant region, decreasing lengths (L—i) /2 of the two connecting regions 
at the ends and decreasing height It appears that at fixed L all these solutions have 
a larger V than the one of the middle type. 
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Figure 5: The functions q{z) in Eq. ([5T|l and L{z^,) in Eg. ([50]) for tlic high temperature metric with 
T — 400MeV using the deformation ((44|) (upper two panels). The lower two panels show the same for 
the low temperature metric with T = 0, now plotted vs z in units of 1/GeV. 



Eqs. (j49p and (j54p diverge for z ^ and have to be renormalized. There are two possibil- 
ities: 

1. Neglect the 1/e terms in (09]) and (|54p and put e = in the remainder; 

2. Subtract the two and put e = in the remainder, i.e., redefine V as V — VqQ- 

We use the first one and the final potential is given by combining Eqs. (fSOl) . (091) and ([Ml) 
so that the equilibrium state always corresponds to the smaller one of V and Vqq . 

Consider first the low T metric, using the deformation ()45p and parameters obtained earlier. 
We plot q{z) and L(2:*) with z and in units of 1/GeV in Fig. [5l The structure in L at 
small is significant and leads to up to three values of V corresponding to one L as seen in 
the plot of the (T-independent) potential V{L) vs L in Fig. [6j Tracing through the curve 
L = L{z^) in Fig. [5] one can trace through the curve V = V{L) in Fig. [HI The almost 
constant branch of V{L) arises from the 2* > 6GeV~^ part of L{z^,). 

We interpret the resulting T = potential as showing a confining part for L < 1.2 fm and 
a constant part for L > 1.2 fm. 

For L < 1.2 fm we can express the potential in Cornell form 

4 ry 

V{L) = --^ + aL. (57) 
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Figure 6: V{L) vs L as computed from eq. ((54)) for the low temperature metric (left, T-independent) 
and F{L, T) vs Lfor the high temperature metric (right) at T = 500, 300, = 186, 100 MeV. Note 
the difference in L-scales. The almost horizontal part of the low T potential comes from the 2;* > 6 
part of the L(z,)-curve in FiglH the horizontal parts of the high T potential are L-independent values 
of Vqq from dUl). 

If we use as ~ 0.36, then we get the AdS radius C? jii^oi') « 0.67 and a 0.168GeV^ or 
\fa w 410MeV. This is in reasonable agreement with recent lattice data giving y/a ~ 460 
MeVga]. Note that T^j^ « 0.46 in this case. 

Beyond 1.2 fm the potential is approximately constant and grows slowly from ~ 0.9 GeV to 
Vqq = 0.961 GeV, obtained by evaluating Eq.([39]) with zq = 00. One can speculate that the 
constant potential is an indication of QCD string breaking of dynamical quarks [17], The con- 
stant part starts at 1.2 fm and its energy scale is about 1 GeV. This should be in principle the 
twice the mass of the lightest hadron. Current lattice data indicates that the string breaking 
occurs about 1.2 — 1.4fm at the energy scale of 1.0 — 1.2 GeV [48j. This is consistent with our 
model calculation. The reason why we see the string breaking is not obvious at first thought 
because it happens only when dynamical quarks are included. However, from the considera- 
tion of linear Regge behavior of meson mass spectrum and the boundary thermodynamics we 
have been lead to a specific form of the function 'w{z) containing a scale related to the mass 
scale of the fundamental matter in QCD. Thus in our model ^ we are already implicitly 
including dynamical quarks from the beginning. In AdS/CFT, light dynamical quarks can 
be added by introducing a spacetime-filling D7 brane [29]. It would be very interesting to 
see in supergravity whether the original metric ansatz ([2]) can be related to this flavor brane 
background. 

Now that we have determined the AdS radius £^/(7ra') ~ 0.67, we can draw the QQ 
potential at high temperatures, again using the deformation (j44p and parameters determined 
earlier. The quantities q{z) and L{z^) with z and z^ in units of zq are plotted in Fig. [5] and 
V{L) in Fig. [6]for T = 400 MeV. It is again useful to trace through L{z^) and see how it maps 
to V{L): one starts from L = 0, F = — cxd, then L and V grow to a maximum value after 
which both again decrease but so that V is almost constant. Finally again L = at = 
at which point V{L = 0) = Vqq. This is L-independent and is shown as the horizontal curves 
in Figini 
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One can observe the following for the high temperature phase: 

1. The overall potential is the lowest of the curves/line and contains the usual conformally 
invariant V ^ —1/L at small L. At some L this reaches Vqq and from this L onwards 
the preferred state is the QQ state with a constant L-independent potential. This is 
the structure observed already in [55] for an undeformed metric. 

2. At very large T, small zq, the large L value of the potential (free energy) Vqq from 
Eg. ([39]) approaches C'^ /{iTa'){—7rT). Although one cannot strictly separate internal 
energy E and entropy S in the free energy of the QQ state, this is compatible with 
E = 0, S = C? I(tioI^ ■ IT « 2.1. Notice that AS ~ 0(1) for simple quark and antiquark 
system. 

3. As the temperature decreases below Tc, a first order phase transition occurs as we have 
seen in sections 3 and 5. One can continue plotting V{L) in this metastable phase, as 
done in FiglU and obtain a linear confining potential. However, the stable phase is the 
one corresponding to the metric ^ which should be used for T <Tc. 

7 Discussion 

In this paper, we have considered a 5d gravity action which presumably could descend from 
a higher dimensional fundamental theory. Following [23j, we proposed a solution of the 
equations of motion that is dual to a QCD-like gauge theory on the 4d boundary. Using 
standard gravity / gauge theory duality correspondence we computed the 4d boundary energy- 
momentum tensor. The gauge dynamics on the boundary is that of a strongly coupled 
fluid modified by a vacuum energy. The bulk geometry was further constrained so that 
the 4d boundary theory shows a first order phase transition. With constraints from the 
hadron spectrum, the parameters of the deformation were fixed. Then we considered the 5d 
bulk transition following the recipe of Hawking and Page [24| and found that the geometric 
transition in the 5d bulk precisely matched with the phase transition on the 4d boundary. 
This is a remarkable correspondence and shows the power of gravity/gauge theory duality 
and the reasonableness of the model we have set up for the dual construction of QCD-like 
theories. The computation of quark- antiquark free energy (potential at T = 0) showed that 
at T = potential contains a confining linear part but also string breaking, indicating that 
the deformation contains effects of dynamical quarks. At T > Tc we saw a deconfined state 
with entropy = 2.1 /pair. 

In general terms, the whole setup started from three unknown functions, the deformations 
of the high T and low T metrics and the dilaton. We chose simple ansatze for the high T 
deformation, h{z) = —cz'^, and for the dilaton, ^{z) = (pz"^. The low T deformation was then 
constrained in various ways and one ended up with the function w{z) plotted in FiglH There 
is some freedom to change this function in the intermediate z range, but not much without 
destroying the coincidence of the 4d boundary and 5d bulk transitions. 

Many QCD-like features treated in this paper are basically IR dominant phenomena, but 
since from gravity/gauge theory duality, the boundary data is critically dependent on the 
boundary behavior of the bulk metric, the metric form in the UV region near the boundary 
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becomes as important. If we solve the equations of motion of supergravity in top-down 
approach, this UV/IR behavior is naturally encoded in the solution, but when we approach 
phenomenologically we have to put this UV/IR behavior by hand. What we saw here is that 
the form of the metric is stringently constrained by this. 

There are several directions of future study. First is to check the our ansatz using various 
physical quantities in QCD and refine the form of w{z). Also it would be highly desirable to 
reproduce the solution in supergravity. Even though this is very difficult task, we have some 
better guiding principle at least. 
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